Noise-resistant control for a spin qubit array 
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We develop a systematic method ol performing corrected gate operations on an array of exchange- 
coupled singlet-triplet qubits in the presence of both fluctuating nuclear Overhauser field gradients 
and charge noise. The single-qubit control sequences we present have a simple form, are relatively 
short, and form the building blocks of a corrected CNOT gate when also implemented on the inter- 
qubit exchange link. This is a key step towards enabling large-scale quantum computation in a 
semiconductor-based architecture by facilitating error reduction below the quantum error correction 
threshold for both single-qubit and multi-qubit gate operations. 



The prospective scalability of spin qubits in semicon- 
ductor quantum dots, along with their demonstrated 
long coherence times and rapid gate operations, make 
them a leading candidate for quantum computing archi- 
tectures. Singlet-triplet qubits [TJ |2J, where the quan- 
tum information is encoded in the zero-projection spin 
states of two electrons in a double quantum dot, are 
particularly promising due to their insensitivity to stray 
magnetic fields and their purely electrical controllability. 
However, precise experimental manipulation of singlet- 
triplet qubits is hindered by two sources of noise in- 
variably present in all laboratory systems: fluctuations 
in the background nuclear spin bath due to long-range 
hyperfine-mediated flip-flop processes (Overhauser noise) 
[3HS], and fluctuations in the electrostatic quantum dot 
confinement potential due to background electrons hop- 
ping on and off nearby impurity sites (charge noise) 0(7]. 
These fluctuations are slow compared to typical qubit ro- 
tation times, and this highly non-Markovian characteris- 
tic can be exploited to suppress their effects by spin echo 
or related dynamical decoupling protocols when one does 
not aim to rotate but instead preserve the qubit state, 
i.e. quantum memory [8HTU]. The ability to achieve very 
long quantum memory times in this way is one of the 
great advantages of semiconductor spin qubit architec- 
tures. The problem is that the dynamical decoupling 
scheme to preserve coherence does not work during the 
qubit gate operations. 

Similar protocols for robust qubit rotations are, there- 
fore, highly desirable so that quantum coherence is pre- 
served during the gate operations. However, singlet- 
triplet qubit control is subject to severe physical con- 
straints that make this task quite daunting. A large 
inter-dot Overhauser field gradient introduced by nuclear 
spin pumping [11H13) (or an actual magnetic field gradi- 
ent due to a proximal micromagnet |14j ) performs ro- 
tations about the x axis of the Bloch sphere, but this 
gradient is not tunable on the time scale of an operation. 
Meanwhile, electrical control of the inter-dot exchange 
coupling by tilting the double-well potential [T] leads to 



rotation about the z axis of the Bloch sphere, but the 
sign of the coupling is fixed. The only rapidly tunable 
element is the magnitude of the exchange coupling, and 
its range is limited by either the singlet-triplet splitting 
of two electrons on a single dot, or, more restrictively, by 
keeping the double-well near a symmetric configuration 
to avoid converting the spin qubit into a fragile charge 
qubit. Thus, one is restricted to positive rotations about 
a limited range of axes somewhere between +x and +£, 
and pulses approximating a delta function are not avail- 
able. This unique set of constraints prohibits straight- 
forward application of control techniques from the NMR 
literature (e.g., [TS1 IS])- The high-fidelity gate opera- 
tions crucial to scalable quantum computation requires a 
totally new approach to the quantum control problem. 

Previously, we have shown that there exists a new 
form of control sequence that respects these constraints 
and eliminates the leading-order single-qubit error due to 
Overhauser noise |17] . Recent numerical work has shown 
that both relevant types of error can be simultaneously 
addressed [T5]. Despite this progress, there remains a 
need for a protocol that corrects both Overhauser and 
charge noise errors while being sufficiently simple and 
flexible for incorporation into multi-qubit operations. 

In this work, we present a method of pulse design that 
achieves this goal, systematically eliminating both errors 
to leading order for any quantum circuit. In the language 
of the NMR community, the task is similar to correcting 
arbitrary quantum gate operations for both amplitude 
and detuning errors simultaneously. This task has only 
recently begun to be addressed in NMR [THl HO] > and it 
is remarkable that the far more restricted case of singlet- 
triplet qubits permits a simple solution, as revealed by 
our method. Furthermore, our new pulse sequences are 
an order of magnitude faster than earlier, less-powerful 
sequences [T7] . We begin by showing how to perform uni- 
versal, robust, single-qubit gates. We then demonstrate 
how to use these single-qubit pulse sequences to generate 
a CNOT gate that possesses the same resilience against 
errors. Finally, we show how to combine the prior two 
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results to generate universal, multi-qubit, dynamically 
corrected operations on a large-scale quantum register. 
Thus, our work forms a complete prescription for com- 
pensating low-frequency noise in ongoing singlet-triplet 
quantum computation experiments. 

The model Hamiltonian within the logical subspace of 
a singlet-triplet qubit is written in terms of the Pauli 
operators a as 



h J (6(f)) 

H {t) = -jVx -\ ^ az 



(1) 



with h — gfiB^-Bz the energy associated with the aver- 
age magnetic field gradient across the double-dot and J 
the positive, bounded exchange coupling. The exchange 
is a function of the energy difference between balanced 
and imbalanced singlet charge states, e, which can be 
controlled dynamically pQ. Since both Overhauser and 
charge noise are typically several orders of magnitude 
slower than gate times, the resulting perturbations about 
h and e (i), Sh and Se, respectively, are treated as random 
constants. 

A single-qubit rotation of angle about axis hi + Jz, 
R(hx + Jz, </>), naively performed by holding the ex- 
change constant over some time interval, results in errors 
Aj (see Supplement for the explicit formulas 21 ), 



U (J, 4>) = exp 
'h 



exp 




J 

2 2 V v^2 + J2 



Here h and Sh are assumed to be independent of J, and 
5 J = 5e dJ £' \j! e )=j arises from fluctuations in the back- 
ground impurity potential and hence in detuning, Se. 

Our general strategy is to construct an identity oper- 
ation such that the error in its implementation exactly 
cancels the leading order error in the original rotation. 
For example, a 2n rotation interrupted by a 2ir rotation 
about a different axis (the SUPCODE identity of Ref. [T7] ) 
allows three degrees of freedom with which to tune the 
error: the two axes and the point of interruption. Here 
we use a similar concept in a more general form to com- 
pensate for all error sources. First, note that a 2m„7r 
rotation interrupted by a general zeroth order identity 
gives a new zeroth order identity with different first or- 
der errors, 

U (jn,rn n n + 6„) ^1 - S^^a^j U (j n , m n ir - 6 n ) 

= J-i£*< n V, (3) 
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This defines a recursion relation for the error of a 
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FIG. 1: (a) Example pulse sequence for it/2 rotation about 
x + z. (b) Parameters vs angle for rotations about x + S. 



"level-rt" parameterized identity, 
U (j n ,m n Tr + n ) ...U (ji,miir + 6i)U (j ,2m 7r) 

X U (jl, TOl7T — 6\) ...U (j 

)• (4) 

The recursion equations are straightforward to generate 
[2"T] . but their explicit form is lengthy and unnecessary 
for the discussion here. 

We find that a level-5 identity contains enough flexi- 
bility to obtain errors which exactly cancel those of the 
naive rotation U (J,(f>) to leading order in both Sh and 
Se. This is the central new result of our work, and we 
choose a simple form for this compensating identity, 



U J,7T 



U (j4, tt) U (j 3 ,7r) U (0, tt) U (ix, tt) U (j , 4tt) 



x U (h , tt) U (0, tt) U (j 3 , tt) U 4 , tt) U [J, tt + 



exp 



J 



Vh 2 + J 2 



-O 



(Sh + Se)' 



(5) 
12 so 



Here we have taken j'5 = J, ^1.2,3,4 = 0, and #5 
that the overall sequence is symmetric and A y +Sy = 0. 
We have arbitrarily set j 2 = for simplicity. The remain- 
ing parameters Jo,Ji,j3j ji are determined by numerical 
solution of the four coupled, nonlinear equations that set 

(5) (5) 

the coefficients of Sh and Se in both A X +5 X and A z +S z 
to be zero |21) . 

An example of this sequence, along with the depen- 
dence of the parameters on the rotation angle, is shown 
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FIG. 2: Gate infidelity of naive (solid) and corrected (broken) 
7r/2 rotations about x + z vs (a) Overhauser field gradient 
fluctuations, Sh/h, and (b) detuning fluctuations, Se/eo. 

in Fig. [T] While there are multiple possible solutions, 
we have constrained our numerical solution to give only 
positive values of the exchange coupling within a typi- 
cally experimentally accessible region. In all cases, our 
restrictions to physically practical sequences still permit 
error compensation. The suppression of error is evident 
in Fig. [2] The corrected sequence infidelity scales as the 
fourth (instead of second) power of the fluctuations. For 
fluctuations on the order of a few percent, our pulse se- 
quence suppresses gate error by two orders of magnitude. 

The results shown here are for the usual empirical 
model at negative e, J (e) oc Jq + Ji exp (e/eo) with J ~ 
so SJ/J = 5e/eo [221 E3] - However, our method does not 
require this assumption, and is equally valid for other 
models of J (e), as we have also verified explicitly, includ- 
ing models where the exchange cannot be tuned all the 
way to zero. The important requirement is simply that 
there is a model so that the correlations between the 8 Js 
at different values of J are known. Also, although we 
have taken the pulse sequence to be piecewise constant, 
we have checked that including a finite rise time only in- 
troduces a perturbation about the parameter values plot- 
ted. In an experimental context, then, an optimization 
of the actual parameters around the ideal pulse described 
here should quickly converge to the adjusted sequence for 
that particular real setup. 

The particular realization of the identity above is not a 
unique solution - we had more parameters than equations 
and simply chose values for which solutions were easy to 
find. Further optimization may yield an even shorter 
solution. However, this identity uses only 147r of total 
rotation and is already quite efficient. (Note that an 
x + z rotation with the original SUPCODE sequence [17] . 
only correcting Sh, naively required 1207T of total rotation 
in the identities, making the current protocol more than 
an order of magnitude more efficient than the original.) 



So far, we have only discussed correcting rotations 
about axes of the form hx + Jz. Given the experimen- 
tal constraints on h and J, this only allows direct access 
to a segment of the positive x-z quarter-plane. Clearly, 
one can build an arbitrary corrected single-qubit rotation 
from a string of the corrected rotations discussed above. 
More efficiently, though, one could also build it from a 
string of uncorrected rotations about the accessible axes 
and then correct the total error with a single compensat- 
ing identity rather than correcting each segment individ- 
ually. We have found that more general rotations can be 
corrected using a level- 6 identity |21) . 

Now we turn to the construction of corrected two-qubit 
gates on a pair of singlet-triplet qubits. Uncorrected two- 
qubit gates have been demonstrated experimentally via 
capacitive coupling [22] . and proposed theoretically via 
exchange coupling [33] [S3] . We will consider the latter 
case with neighboring qubits A, consisting of dots 1 and 
2, and B, consisting of dots 3 and 4. The qubits in- 
teract via an exchange link between dots 2 and 3. In 
recent work, Ising gates are constructed in the absence 
of noise such that phases due to static local Overhauser 
fields cancel, and only a state-dependent phase from the 
exchange pulse survives [24], 

U xx (a) = % tt) (z, tt) C 23 (I) 
xRW (z,tc)RW (^TrJCaag) 

= exp(i~*WaW)+0(5h,5e), (6) 

where R^ A / B ^ (f,9) denotes a rotation of qubit A/B by 
6 about f on the singlet-triplet Bloch sphere. C23 (a/2) 
denotes the application of a pulse to the inter-qubit link 
such that in an S z = subspace it would act as a 2ir rota- 
tion about some axis. This is required to avoid swapping 
anti-aligned spins on dots 2 and 3, which would cause 
leakage out of the logical qubit subspace. The choice of 
axis is determined by fixing the total inter-qubit pulse 
area, J dt J23 (t) — a/2, which causes the desired relative 
phase to be acquired by two-qubit states with aligned 
electron spins on dots 2 and 3. 

By using our compensated pulse sequences to perform 
the R( a / b } (z,tt) operations - and performing them on 
A and B with equal durations such that both qubits are 
idle for the same amount of time - single-qubit errors 
on the rhs of Eq. ^ are eliminated to first order. For 
equal gradients on both qubits, the equal time condition 
is trivially satisfied. More generally, the tt rotations on A 
and B would have different durations and one would pad 
them with corrected 2-7T rotations about different axes 
(i.e., having different durations) to compensate [21] . 

Also implementing the C23 (a/2) with our pulse se- 
quences suppresses leakage error, but the angle a is still 
sensitive to charge noise induced exchange amplitude er- 
ror. This is because our sequences only correct rotations 
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FIG. 3: Quantum circuits for (a) Overhauser-corrected and (b) fully-corrected two-qubit Ising gates. 



within the S z — subspace of the SU(2)® space on 
which C23 acts. The problem at this stage, though, is 
one familiar from NMR contexts - pulse length error in 
Ising-coupled spins [25]. Defining a tilted version of the 
Ising gate, 

U' xx (a, 4>) = R {A) (z, ~<t>) U xx (a) (£, 0) , (7) 

the inter-qubit exchange amplitude error can be cor- 
rected using a BB1 sequence [13125]. Typically we only 
report angles modulo 2n, but because the error is pro- 
portional to a and using our composite pulses results in a 
large angle a — Sir + mod (a, 2tt), we must use a slightly 
modified version of the BB1 sequence: 



U' xx [6 u <h) U' xx (20 1; 30!) U' xx (ft, 0i ) U xx (a) 



exp 



(.a 
V 2 



(A) a (B) 



O 



(Sh + SeY 



(8) 



where Q\ = 8n + n, <j>i — arccos (—a/4 (8ir + n)) [27], and 
all rotations are corrected composites. For a — 8ir + n/2, 
this gate is equivalent to a CNOT gate up to single-qubit 
operations [23][2S]. In practice, we also perform corrected 
identities on qubit B during the rotations on A in Eq. Q 
to avoid reintroducing single-qubit errors. 

We sketch the implementation of the Ising ir/2 gate 
in Fig. [3] We have displayed the case of a linear gra- 
dient, after trivial contractions of sequential single-qubit 
operations in Eq. Q. Each line represents a quantum 
dot, which may be linked to its neighbor by an exchange 
pulse. A singlet-triplet qubit then is denoted by a pair of 
lines which are understood to have total spin projection 
zero. The error suppression for a CNOT gate is shown 
in Fig. [4] For the sake of display, we have reduced the 
number of parameters by assuming a uniformly fluctuat- 
ing linear field gradient across the four dots, as well as 
identical detuning fluctuations on each link. 

As is always the case with composite pulse schemes, 
the cost of the correction is longer gate time. Currently, 
our shortest implementation of the corrected CNOT gate 
requires 20 composite pulses, corresponding to ~ 3007T of 
total rotation. While certainly challenging, this is well 
within the realm of possibility given that our sequence 
suppresses decoherence during its operation. Recall that 
statistical fluctuations of the control Hamiltonian lead to 
rapid decoherence on the free induction time scale, T 2 . 
Dynamical decoupling can extend coherence to a much 



longer time T2, though, where dynamical fluctuations be- 
come important. Roughly speaking, gate error scales like 
the ratio of gate time to coherence time. Since typically 
T 2 > 10 4 T 2 * in semiconductor spin systems [Tl 121151 [25]. 
it is well worth using longer gate implementations that 
simultaneously perform dynamical decoupling such that 
the relevant coherence time is close to T 2 rather than T 2 . 

Since each experiment will have its own preferred re- 
gion of control space, further optimization would be pre- 
mature, but future work will almost certainly reduce the 
length in any case. Also, it may well be possible to con- 
struct a "one-shot" correction for the CNOT gate, as we 
have done for the composite single-qubit rotations, rather 
than correcting at each intermediate stage. These exten- 
sions, however, are beyond the scope of the present work. 

Finally, we discuss application to an array of exchange- 
coupled qubits. Our method facilitates operations below 
the quantum error correction threshold, a vital require- 
ment for scaling up. Arbitrary corrected multi-qubit cir- 
cuits can be implemented similarly to Fig. [3] as long as 
there exist single-qubit corrected identities of sufficiently 
variable duration to protect the many idle qubits during 
the application of nontrivial gates to other qubits. We 
find that level-6 identities can be used to cover the entire 
range of relevant idling times [2T] . An additional benefit 
of using our sequences in large systems is that errors due 
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FIG. 4: Gate infidelity [5T] of uncorrected (solid) and cor- 
rected (dashed) CNOT operations vs (a) Overhauser field gra- 
dient fluctuations, 8h/h, and (b) detuning fluctuations, Se/tQ. 
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to small spatial inhomogeneities in the control Hamilto- 
nian ([I]) are also suppressed. 

In summary, we have shown above that the relevant 
types of error in experimental singlet-triplet spin qubit 
manipulation can be eliminated to leading order by a new 
composite pulse sequence. We show how to apply this se- 
quence for single-qubit, two-qubit, and multi-qubit sys- 
tems, opening a path towards scalable, universal quan- 
tum computation in a noisy solid-state environment. 

This work is supported by IARPA and LPS. 
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Supplementary material 



I. ERROR TERMS IN U(J, </>) AND UNCORRECTED IDENTITY 

In this section we present the explicit form of the error terms in U(J, (/>), which is defined in Eq. 
We start with the error terms in Eq. (|2|), which, to the first order of Sh and 5 J, are 

xu h 2 4> + J 2 sin0 foJ(0-sin</>) 
A x = dh ^-pr+dJ — -, (b-la) 

2(h 2 + J 2 ) 3/2 

, , -h(l — cos 6) ,_ „_ , 

y = Sk - ™ + 6J 2 W + JV (S_lb) 



2(/i 2 H 


- J 2 ) 3 / 2 


J (cos 


0-1) , 


1 2 (h 2 ■ 


f J 2 ) 


hJ(<f>- 


— sin </>) 


2 (ft 2 4 


- .Pf' 2 



(./ J o • // J >ino) 

A z = Sh ^ ^ +S J^ -%L. S-lc 

• ■■' 2(ft 2 + J 2 ) 3/2 

We aim to cancel these first-order error by supplementing the rotation with a level-n uncorrected identity /, 

T(-)= I -iY^5^a i , (S-2) 

i 

such that Aj + #j = for i — x, y, z. 



() 



The uncorrected identity we use in our work is denned recursively in Eqs. ^ and Q. It is therefore useful to 
establish the recursion relation for the coefficients. At zeroth order, is simply U (jo, 2mo7r), therefore 

c( ) _ ra 7rfe (hSh + jpSjo) 

4 0) = 0, (S-3b) 
s ( ) _ "-Wo {hSh + j Q Sjo) 

According to Eqs. Q and Q, given j„, m n and <5 l -™' ) can be generated recursively as 

S ( n ) = m n wh(hSh + j n Sj n ) ^ (n _ x) ft 2 + (-1)"'" j 2 cos6>» ^„_ 1} ,_ 1 s,„„ j n sin 6>» 

(/> 2 +j 2 ) 3/2 x h 2 + J 2 y v/^+jI 

+ ^- 1) fejn(l-(-ir-COSfl„) 

e) = (-l) m » (s£-vifih= + 8f-V cos e n + 5<?-V h8in6n ) , (S-4b) 



5 (n) = m n irj n (h5h + j n 5j n ) s{n _ 1} hj n (l-(-l) m " cos9 n ) _ s{n _ 1} /_^m„ 



/i sin 8„ 



-i) m - h 2 

h 2 +.] 2 



+ f (,-l) in + (-l) C 0S g„ _ (s _ 4c) 



For convenience, from here onward we take ft. = 1 as our energy unit. We also denote g(J) — \ j(e)=J so that 
6 J = g{J)Se. 



II. ROTATION AROUND AXIS x + Jz 



This section discusses how one can correct a rotation of angle (j) around an axis determined by x + Jz, U(J, <f>). 
According to the main text [see Eq. J5J], our target is to design an uncorrected identity /, satisfying 



U (J, 4)1 = e lx R(x + Jz, <j>)ll + (Sh + Sef } 



(S-5) 



where x is an unimportant phase factor. 

We correct this rotation using a level-5 identity. We have also chosen the pulse to be symmetric, namely 6*5 = — </>/2 
and #1,2,3,4 = 0. In this case, the a y term would automatically vanish in the final rotation (see the "Methods" section 
of Ref. |17j). and we only need to solve four coupled nonlinear equations. 

We can absorb the outmost level into U (J,<p) and are left with a level-4 identity whose parameters are fixed by the 
conditions 



<5l 4 > = - 



27r + 0-2J 2 sin| J 5 (J)(27r + + 2sinf) 
A 5h — 



5^ = - 



2(1 + J 2 ) 3/2 
J ( 2tt + 6 + 2 sin 



Sh- 



2 (1 + J 2 ) 3/2 
g(J) \j 2 (2tt + 4>) 



5J, 



2 sin 



SJ. 



2(1 + J 2 ) 3/2 2(1 + J 2 ) 3/2 

As explained in the main text, we can construct the identity as 

1^ = U (j 4 , ?r) U (j 3 , tt) U (j 2 , tt) U (ji,7r) U (jo, 4tt) U (ju tt) (j 2 , tt) U (j 8 ,7r) (7 (j 4 , tt) 



(S-6a) 
(S-6b) 

(S-7) 



Since we only have four equations, we only need four parameters to determine. Here we choose j'2 = 0- (We have also 
verified that our method does not require this assumption to work, namely one can choose a different value of ji , or 
one can instead fix another variable, e.g. 
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Supplementary Figure SI: Parameters v.s. angle for rotations about z. g(J) = J/eo- The shaded area around 0.6-7T < (f> < 0.97T 
indicates the range where j'3 becomes negative, rendering the solution unphysical. 



III. ROTATION AROUND i-AXIS 



The zeroth order z-axis rotation, R(z, <j)), is achieved by [17j : 

R(x + z, -k)R(x, <fi)R(x + z, 7r). 



(S- 



At first glance, since each of the three terms in Eq. ( S-8 ) can be implemented using the pulse of Sec. [TTJ it is 
straightforward to implement a dynamically corrected z-rotation. However, this would require 40ir ~ 50tt of sweeps 
around the Bloch sphere, we are interested in further optimizing it. We shall try to do a "one-shot" correction, namely 



correcting the pulses of Eq. (S-8) with only one identity. 



One trick to be employed here is that we shall insert the identity between R(x, <p) and R(x + z, tt) but not at the 

(S-9) 



right end of Eq. ( S-8 1 , so the corrected pulse looks like 

U (J = 1, n)U(J = 0, cp)I (6) U(J 



1,7T 



As in Sec. [nj we absorb the outmost level of that identity into U( J — 0, </>) so that the corrected pulse is 

U(J = l,n)U ( J = 0,tt+ I ) I (5) U ( J = 0,tt+ i ) U(J = 1,tt). 



(S-10) 



Therefore if we implement 1^ in a symmetric way then the entire sequence is also symmetric, assuring us that the 
<J V component will automatically vanish and we are only left with four rather than six equations to solve. This is why 
we place the uncorrected identity in the middle of the sequence. In fact, an uncorrected identity operation can be 
placed anywhere and one may simply choose a location which is most convenient. 

We will not explicitly carry out the algebra here, but we find the following sequence already accomplishes our 
purpose: 



U(J=l,7r)U(j 5 = 0,2tt+ |)C/(j4,7r)C/(j 3 ,^)[/(j 2 ,7r)C/(j 1 = 0, n)U{jo, 4tt) 
x U{h = 0,ir)U(j 2 ,ir)U(j 3 ,ir)U(j 4 ,Tr)U(j 5 = 0, 2tt + ^)U(J = 1,tt) 



(S-ll) 



The parameters jo. 2, 3, 4 are given in Supplementary Fig. SI Note that for a small range of </> (0.6-7T ~ 0.97r), j 3 becomes 



negative, making the solution unphysical. Although this problem can be easily overcome by duplicating pulses that we 
can physically do, it is more natural to stick to a "one-shot" correction to keep it optimal. In fact, we have explicitly 
verified that, setting j'2 = 3a = (instead of ji = js = here) while solving for jo, ji, J3, 35 will generate a physical 
solution to the pulse sequence covering this <p range. 

We now consider the total rotation angle of the pulse sequence Eq. ( S-ll ). We see that it requires around 187T ~ 20ir 



of rotation on the Bloch sphere, and is more than a factor of two shorter than the naive speculation at the beginning 
of this section. 
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Supplementary Figure S2: Parameters v.s. angle for rotations about y. g(J) = J/eq. 



IV. ARBITRARY ROTATION 



It is well-known that an arbitrary rotation can be decomposed as 

R(x, 4>a)R{z, 4>b)R{x, C ). 
As we already know from Sec. |III[ this decomposition can be written as 

R(x, 4> a )R(x + z, n)R(x, (f>b)R(x + z, tt)R{x, (j> c ) 



(S-12) 



(S-13) 



up to a phase factor. Again, although we have the ability to correct each term in Eq. (S-13), we would prefer a 
"one-shot" correction at the cost of introducing a slightly higher level of uncorrected identity operation. 

We note that for an arbitrary rotation, a y terms are in general present, so typically we have to solve all six equations. 
The identity is inserted to the sequence as 



U( J = 0, <j> a )U{ J = 1, n)IU( J = 0, <j> b )U( J = 1, ir)U( J = 0, </> c ). 



(S-14) 



We still insert / in the middle of the sequence with the hope that when the rotation axis is on the x-z plane but is 
not implementable using methods described in Sees. |n] and |III| (for example in the second quadrant of the plane, or 
in the first quadrant but very close to the z-axis), this entire pulse sequence can still be made symmetric so one is 
to solve four instead of six equations. In the work presented here we use the sequence in Eq. (S-14) since we already 



find physical solutions for a vast range of rotations. However one should note that this is not necessary, and one can 
seat the uncorrected identity (or identities) at any position in the sequence as long as it provides sufficient degrees of 
freedom to find physical solution to the array of six coupled nonlinear equations. 

We illustrate our sequence with a R(y, 4>) rotation, which is an important operation used, for example, in converting 
the two-qubit Ising gate to a CNOT gate [52] [55] . Obviously this would have a y terms in the final sequence and one 
need to solve six equations. We therefore consider a level-6 identity with 9q and jo through j§ to be determined. (We 
take j2 = to keep the number of unknown variables six.) The pulse sequence for a corrected y-rotation is 

U(J = 0, a = —)U(J - 1, n)U(j 6 , TT - e 6 )U(j 5 ,w)U(j 4 , n)U(j 3 ,w)U(j 2 = 0, ir)U{h^)U{j 0: 4ir) 

2 (S-15) 
x U(j U 7r)U(j 2 = O,7r)U(h,7r)U(u,n)U(j 5 ,Tr)U(j 6 ,7r + 6 )U(J = O,0 fc = <j>)U(J = l,n)U(J = O,0 C = ~), 

and numerical results of jk and 9 e for a range of rotation angles are shown in Supplementary Fig. |S2[ We remark 
that the pulse sequence Eq. (S-15 1 sweeps a total angle of 20ir ~ 22ir around the Bloch sphere. Comparing to a naive 
correction of an x-z-x sequence which would cost 40-7T ~ 50-7T, this is again a factor of two improvement. 



V. IDENTITY OPERATION 

The identity operation /, sometimes called the NOOP operation, also plays a role in applying our method to two- 
qubit and multi-qubit system. First, at the two-qubit level it provides a fundamental ingredient of the entangled 



9 



gate, defined as C23 (f ) [cf. Eq. (|6| and Fig. |3|a)]. Here a is related to f dtJ 2 3 (t), which, in turn, determines the 
argument of the Ising gate U xx (a). Second, since the Overhauser field gradient is always present during the entire 
duration of performing a quantum algorithm, if only a subset of the multi-qubit system is being gated, the remaining 
qubits must carry out corrected identity operations in order to reduce vulnerability to noise. In this case, the timing 
is crucial. For example if one of the qubits is being rotated with time duration T, then the rest of the multi-qubit 
system must perform corrected identity operation with the same time duration T [cf. Fig.[3^b)]. On the other hand, 
if one qubit is being gated with time t\ while another t 2 , then they must be supplemented by identities with time 
duration T — t\ and T — ti respectively while others perform identity with time T. 

These considerations mean that we must find identity operations which can generate a range of time duration as 
well as J dtJ it). We managed to do this with the following pulse sequence: 

U{J,2n)U{jB,ir)UUi = 0,ir)Ufo,*)U{fr = Q^)U{h,it)U[j ,^) ( . 

xU(juir)U{j2 = Q,n)U(j3,*)Ufa = 0,«)U(js,ir)U{J,2n) 



The pulse sequence of Eq. (S-16) generates identity operations with T between 11/h and 2A/h by varying J, 
from which identities with any desired time T > 11/h can be derived. We have used these identity operations in 
constructing the two-qubit gates schematically shown in Fig. [3](b). This pulse sequence at the same time generates 
the integrated value of J(t), J dtJ(t), between 10 and 20 (note that this corresponds to 20 < a < 40). Therefore 
in Fig.||b), U xx (17n/2) and U xx (9ir) can be generated using this sequence. U xx (18w) can either be generated by 
duplicating U xx (9ir), or, more optimally, be generated by a sequence with one higher level of identity operation, for 
which we shall not present details here. 

VI. DEFINITION OF FIDELITY 

Figs. [2] and |4] show the performance of our pulse sequences in terms of the infidelity (one minus the fidelity, F). For 
a single qubit the fidelity is straightforwardly defined [Slj. To quantify the fidelity of the two-qubit gate, we adopt the 
result of Ref. [32]. However, note that due to presence of the leakage subspaces in our system, our evolution operator 
is not trace-preserving. Therefore Eq. (32) of Ref. [32] must be extended as 
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^Tr[cr (g) a U f <T ® a Uj] + \ V Tr|Vo> ® a^Ufa^ (g> o v U)\ 

0</i<3 
0<y<3 



(S-17) 



where 00 is the 2x2 identity matrix, and X®Y denotes an operator entirely within the computational subspace that 
acts as X on the first qubit and as Y on the second qubit (for example Co ® <Jo is the projector into the computational 
subspace), V is the desired operation with identity in the leakage subspace, and Uf is the actual time evolution at 
the conclusion of the composite pulse sequence. 
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